Anomalous Dynamic Scaling in Locally-Conserved Coarsening of Fractal Clusters 
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We report two-dimensional phase-field simulations of locally-conserved coarsening dynamics of 
random fractal clusters with fractal dimension D = 1.7 and 1.5. The correlation function, cluster 
perimeter and solute mass are measured as functions of time. Analyzing the correlation function 
dynamics, we identify two different time-dependent length scales that exhibit power laws in time. 
The exponents of these power laws are independent of D; one of them is apparently the "classical" 
exponent 1/3. The solute mass versus time exhibits dynamic scaling with a D-dependent exponent, 
in agreement with a simple scaling theory. 

PACS numbers: 61.43.Hv, 64.75.+g 
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Many non-equilibrium systems develop morphological 
instabilities and ramified growth at an early stage of the 
dynamics, exhibit coarsening at an intermediate stage, 
and finally approach a simple equilibrium. A classic ex- 
ample is diffusion-controlled systems, such as deposition 
of solute from a supersaturated solution and solidification 
from an overcooled melt. The stage of morphological in- 
stability has been under extensive investigation 0-Q]. If 
some noise is present, fractal clusters (FCs) similar to 
diffusion- limited aggregates (DLA) can develop ||. If 
the total amount of mass or heat is finite, the subse- 
quent dynamics are dominated by surface-tension-driven 
relaxation (coarsening) [||-[7| . Coarsening of FCs in sys- 
tems with conserved order parameter, apart from being 
interesting in its own right, exemplifies a more general 
problem of phase ordering: emergence of order from dis- 
order, following a quench from a disordered state into a 
region of phase coexistence ||[|. This example is non- 
trivial because of the long-ranged, power-law correlations 
intrinsic in FCs. The role of long-ranged correlation in 
phase ordering dynamics has been under debate in con- 
nection to dynamic scale invariance (DSI), a major sim- 
plifying factor in theory 0. DSI presumes that there is, 
at late times, a single dynamic length scale l(t) so that 
the correlation function C (r, t) approaches a self-similar 
form g[r/l(t)] |1. In locally-conserved systems (model 
B) l(t) is expected to show dynamic scaling (by which 
we simply mean a power law in time) with "classical" 
dynamic exponent 1/3 

Because of complexity of phase-ordering dynamics, 
DSI has not been proven, except for a very few sim- 
ple models ||. However, there is extensive evidence, 
from experiments and simulations, supporting DSI in 
conserved systems with short-ranged correlations. Re- 
cently, phase-field simulations of coarsening of a globally- 
conserved, interface-controlled system with long-ranged 
correlations were performed [0|. DLAs with fractal di- 
mension D = 1.75 served in Rcf. Illll as the initial con- 



ditions for the minority phase. Notice that a FC is 
a particular case of systems with long-ranged correla- 
tions. D is equal to the exponent a that appears in 
the power-law decay of the correlation function at t = 0: 
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where d is the Euclidian dimension 



C(r, t = 0) 

P). The simulations JTTJ have not found any deviations 
from DSI and yielded the "classical" value of the dynamic 
exponent for the globally-conserved model, which is 1/2. 

The results ]ll| stand in sharp contrast with simulation 
results on locally-conserved coarsening of DLA clusters, 
where breakdown of DSI was observed . Breakdown of 
DSI is related to the fact that the upper cutoff L of the 
FCs remains almost constant in the process of coarsening 
While interpreting this finding, one should distin- 
guish between two regimes: diffusion-controlled: Id < L, 
and Laplacian: Id > L, where Id ~ t 1 / 2 is the diffusion 
length. It is clear that L should remain almost constant 
in the diffusion-controlled regime, as interaction between 
far-lying parts of the cluster is exponentially small in 
this case. The inequality Id < L was satisfied in the 
simulations of Refs. H@], and it is satisfied in the sim- 
ulations reported in the present work. However, even 
in the opposite limit of Laplacian coarsening one can ex- 
pect L to remain almost constant (and breakdown of DSI 
to persist), this time because of Laplacian screening of 
transport. An additional non-trivial aspect of locally- 
conserved coarsening of DLA clusters is the following. 
Though not dynamically scale- invariant, this coarsening 
process was shown to exhibit dynamic scaling (with "un- 
usual" exponent 0.21 - 0.22) ||-[|. This implies hidden 
simplicity of a more complex nature than DSI. 

We report phase-field simulations of locally-conserved 
coarsening of 3 sets of random FCs. The clusters have 
fractal dimensions D — 1.7 and 1.5. Our motivation was 
to check whether the "anomalous" scaling persists for 
clusters different from DLA, and to investigate its pos- 
sible dependence on D. The simulations confirm break- 
down of DSI in all cases. We identify three power laws in 
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time. Two of them correspond to dynamic length scales 
found from the dynamics of C{r, t) at small and inter- 
mediate distances. The first length scale, with exponent 
0.21 — 0.22, is the same as observed earlier. Surpris- 
ingly, it does not show any dependence on D. The second 
length scale has exponent close to 1/3, and we suggest a 
simple interpretation for its appearance. The third dy- 
namic exponent is found in the time dependence of the 
"solute mass", and it is /^-dependent. We suggest a sim- 
ple scaling theory for the second and third exponents. 

Here is a brief description of our simulation techniques. 
We employed a hierarchical algorithm |0| to build ran- 
dom FCs with tunable fractal dimension. The algorithm 
starts with a set of 2™ square particles grouped into pairs. 
These pairs are then grouped into pairs of pairs, and so 
on, in an iterative procedure. After n iterations the final 
aggregate of 2™ particles is obtained. The desired frac- 
tal dimension D is achieved by an appropriate mutual 
position in which the aggregates are sticked together in 
each iteration . The clusters obtained in this way are 
"reinforced", by an addition of peripheral sites as sug- 
gested in Ref. ||, to avoid breakup at an early stage of 
coarsening. The quality of obtained FCs is controlled 
by computing the correlation function (see below) and 
checking the quality of the power law. 

A standard phase-field model for locally-conserved 
coarsening is the Cahn-Hilliard equation 

^+ 1 -V 2 (V 2 u + u-u3)=0. (1) 

This equation was discretized and solved numerically. 
Two different numerical schemes were employed to ad- 
vance the solution in time. The first scheme used an 
explicit Euler integration. The simulated domain f2 was 
a two-dimensional box 1024x1024 with no-flux bound- 
ary conditions (zero normal component of Vm at the 
boundaries). The grid size was Ax = Ay = 1, the time 
step At — 1/25 ensured numerical stability. The time 
range of these simulations was < t < 5000. The sec- 
ond scheme used a semi-implicit Fourier spectral method 
The simulation box was 1792x1792, with periodic 
boundary conditions. The grid size was Ax = Ay = 1 
that corresponds to 1792 Fourier modes in each dimen- 
sion. The time step was At = 0.3, the total time range 
was < t < 15000. Overall, three scries of simulations 
were performed. In two of them, with the explicit scheme, 
10 FCs with D = 1.7 (series A) and 10 FCs with D = 1.5 
(series C) served as the initial conditions for the minority 
phase u = 1. In addition, 7 FCs with D = 1.7 (series B) 
were simulated with the spectral scheme. The minority 
phase area fractions were 10.5% (series A), 7.5% (series 
B) and 5.8% (series C). Notice that the initial cluster 
area in series B was almost twice as big as that in series 
A, while D was the same. 

Eq. (fil) with either no-flux, or periodic boundary con- 



ditions obey conservation laws: Iq = J J it(r,t)dr = 

n 

const and Ii = J J r u(r, t) dr — const. Our numeri- 

o 

cal schemes preserve Jo exactly. We checked that they 
also preserve Ii with a very high accuracy: better than 
5 x 10 -3 % up to t — 5000 (explicit scheme), and better 
than 0.05% up to t — 15,000 (spectral scheme). Ad- 
ditional tests of the two schemes included obtaining the 
stationary kink solution of Eq. ([!]) and observing the 
"classic" scaling with exponent 1/3 for the initial condi- 
tion in the form of "white noise" . 



FIG. 1. Images of locally-conserved coarsening of a FC 
with D = 1.7 at times t — (upper left), 495 (upper right), 
1963 (lower left) and 4920 (lower right). 
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FIG. 2. Same as in Fig. 1, but for D = 1.5. 



The cluster was identified as the locus of B(r,t) = 1, 
where B(r, t) — 1 for it(r, t) > and zero otherwise. Fig- 
ures 1 and 2 show snapshots of the coarsening dynamics 
observed in series A and C. One can see that larger fea- 
tures grow at the expense of smaller ones. The global 
structure of the clusters remains unchanged. In particu- 
lar, the cluster size is almost constant, in a marked con- 
trast to the globally-conserved case [0. A significant 
decrease of L in time is a necessary prerequisite of any 
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fractal coarsening that obeys DSI ]6|, |1 l| , |14| , pL5| , therefore 
its absence clearly implies that DSI is broken. One can 
also notice a more pronounced breakup of the cluster for 
D = 1.5. Still, even in this case the large-scale mass 
distribution in the cluster does not change much. 

To characterize the coarsening dynamics, several quan- 
tities were sampled and averaged over the initial condi- 
tions in each of the three simulation series: 

1. (Circularly averaged) equal-time correlation func- 
tion, normalized at r=0: 



C(r,t) 



< p(r' + r,t)p(r',t) > 
<P 2 (r',t) > 



2. Cluster perimeter P(t) |lq| . 

3. Cluster mass M(t) = J J p (r, t) dr . 

B(r,t) = l 

4. "Solute" mass outside the cluster 
M s (t) = J J p(r,i)dr. 

S(r,t)=0 

Here an auxiliary "density" field p(r,i) = (1/2) [u(r, t) + 
1] is introduced that varies between and 1. Also, Iq = 
const implies that M(t) + M s (t) = const. 
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FIG. 3. Equal-time pair correlation function at different 
times for simulation series B: t = (solid), t = 658.8 (dotted), 
t = 2787 (dashed-dotted) and t = 15000 (dashed line). The 
inset shows a log-log plot of the same function. The dynamic 
length scale h(t) (the knee- like feature) is apparent. 

Figure |^ shows C (r, t) at different times for simulation 
series B (D — 1.7). Coarsening occurs at small and inter- 
mediate distances, while the large-r tail remains almost 
unchanged. The fractal behavior is observed only in the 
"frozen" tail that shrinks in time. On the typical scales 
of coarsening, the cluster is not fractal anymore. This 
implies breakdown of DSI at D < 2. The same type of 
behavior of C(r, t) was obtained for the locally-conserved 
coarsening of DLA clusters W$\. At small distances, 
C(r,t) behaves linearly with r: C(r,t) ~ 1 — r/l\{t), 
as expected from the Porod law ||. This asymptotics 
yields dynamic length scale h(t). One can also see a 
knee-like feature at intermediate distances that develops 



in the course of time. This feature yields dynamic length 
scale li. We define I2 as the maximum value of r for 
which the relative difference between C(r, 0) and C(r,t) 
is not greater than 3% (the results are insensitive to the 
exact value of this threshold). Notice that, by the end 
of the simulation time, the diffusion length Id is already 
greater than 10 2 . Still, we do not see any signature of Id 
in the shape of C(r, t) at large distances. We expect that 
the tail of C(r,t) will remain "frozen" until late times, 
and DSI broken, in the Laplacian regime as well. 

Figures ^ and [| show the measured quantities l\, I2, P 
and M s , averaged over the initial conditions, versus time 
for simulation series B and C, respectively. We fitted 
these time dependencies by power laws with exponents 
a\,/3, a.2 and 7, respectively. According to Porod law 
||, one expects ct\ — ~ct2- All simulation results are 
summarized in Table 1. a\ and 012 are close to the pre- 
viously reported values Remarkably, they do not 
depend on D. The exponent (3 is close to the classical 
exponent 1/3. The appearance of the classical exponent 
in a situation with broken DSI requires an explanation 
(see below). The solute mass exhibits dynamic scaling 
with a Z?-dependent exponent 7. 

Now we report a simple scaling theory for the expo- 
nents [3 and 7. It is based on the fact that, at subdiffu- 
sive distances, r < ld(t), the Cahn-Hilliard dynamics are 
reducible to a sharp-interface model of Laplacian coars- 
ening HQ. The process of Laplacian coarsening can 
be described as follows. Small branches of the cluster 
shrink and disappear, and their material is reabsorbed 
by larger branches. Because of Laplacian screening, reab- 
sorption occurs locally (therefore, the mass distribution 
on large scales remains unchanged). The reabsorption 
events cause undulations of the interfaces of the branches. 
Assuming DSI on length scales h{t) < r < ld(t), we see 
that, by time t, the characteristic undulation lenght is of 
order t 1/3 . This gives a natural explanation to the dy- 
namic length scale I2 ~ t 1 / 3 observed in the simulations. 
The Gibbs-Thomson condition at the interface 0,0 im- 
plies a solute density of order < -1 / 3 in the "contaminated" 
region within a distance Id from the cluster. Taking into 
account the (preserved) fractal structure of the cluster at 
distances large compared to Id, we estimate the contam- 
inated area as A d ~ l\ {L/l d ) D oc <( 2 ~ D )/ 2 . Multiplying 
the "solute density" i -1 / 3 by Ad yields the solute mass: 
M s (t) oc t^ 3D ^ 6 . The last row in the table shows the 
theoretical exponent 7^ = (4— 3D)/6. A good agreement 
with simulations is seen. Notice that jth changes sign at 
D = 4/3. At D < 4/3 reabsorption should effectively 
stop, and the cluster should continue dissolving. We per- 
formed a single simulation with a random FC with D=1.3 
that indeed showed a very slow but persistent increase of 
M s with time. Notice, however, that FC with a smaller D 
also implies, in our simulations, a smaller area fraction 
of the minority phase. At too small area fractions the 
minority phase can dissolve completely, independently of 
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D. Therefore, a more careful investigation is needed in 
order to distinguish between these two effects. 
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FIG. 4. h (a), l 2 (b), P (c) and M s (d) vs. time for 
D = 1.7. The solid lines are power-law fits. 



In summary, the simulations support earlier results 
on breakdown of DSI in locally-conserved coarsen- 
ing of systems with long-ranged correlations and reveal 
new signatures of hidden simplicity in these systems. We 
confirm the value of the anomalous dynamic exponent 
«i ~ 0.21 — 0.22. Surprisingly, it does not show any 
dependence on D. Two additional dynamic exponents 
(3 ~ 1/3 and 7 = 7(D) are identified. They can be inter- 
preted in terms of Laplacian coarsening at subdiffusivc 
distances r < ld(t), combined with DSI at intermediate 
distances h(t) < r < ld{t). 
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FIG. 5. Same as in Fig. |, but for D = 1.5. 



TABLE I. Dynamic exponents found in the simulations. 
Also shown is 7th, theoretical prediction for 7. 



Series 


A 


B 


C 


D 


1.7 


1.7 


1.5 


Ql 


0.21 ±0.01 


0.22 ±0.01 


0.21 ±0.01 


Q2 


-0.21 ±0.01 


-0.21 ±0.01 


-0.21 ±0.01 


P 


0.32 ±0.01 


0.32 ±0.01 


0.30 ±0.01 


7 


-0.16 ±0.005 


-0.16 ±0.005 


-0.09±2- 10" 4 


1th 


-0.18 


-0.18 


-0.09 



The errors are estimated by shifting the time intervals of fitting. 
The errors in each fit are much smaller. 
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